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ABSTRACT 

We present a broad-band timing analysis of the accreting white dwarf system MV 
Lyrae based on data obtained with the Kepler satellite. The observations span 633 
days at a cadence of 58.8 seconds and allow us to probe 4 orders of magnitude in tem- 
poral frequency. The modelling of the observed broad-band noise components is based 
on the superposition of multiple Lorentzian components, similar to the empirical mod- 
elling adopted for X-ray binary systems. We also present the detection of a frequency 
varying Lorentzian component in the lightcurve of MV Lyrae, where the Lorentzian 
characteristic frequency is inversely correlated with the mean source flux. Because 
in the literature similar broad-band noise components have been associated to either 
the viscous or dynamical timescale for different source types (accreting black holes or 
neutron stars), we here systematically explore both scenarios and place constraints on 
the accretion disk structure. In the viscous case we employ the fluctuating accretion 
disk model to infer parameters for the viscosity and disk scale height, and infer un- 
comfortably high parameters to be accommodated by the standard thin disk, whilst 
in the dynamical case we infer a large accretion disk truncation radius of ~ IORwd- 
More importantly however, the phenomcnological properties between the broad-band 
variability observed here and in X-ray binaries and Active Galactic Nuclei are very 
similar, potentially suggesting a common origin for the broad-band variability. 

Key words: accretion discs, stars: binaries: close, stars: novae, cataclysmic variables, 
stars: MV Lyrae, black hole physics, stars: oscillations. 



1 INTRODUCTION 

Compact interacting binaries (CBs) are close binary systems 
usually consisting of a late-type star that transfers material 
onto a black hole (BH), a neutron star (NS) or a white dwarf 
(WD) via Roche-lobe overflow. With an orbital period on 
the order of hours, the donor star transfers material through 
the LI point, which forms an accretion disk surrounding 
the compact object. The dynamics and physics governing 
the flow of matter accreting onto the compact objects is, 
however, still debated. The accretion disks in BH and NS 
binary systems (X-ray binaries or XRBs) emit most of their 
radiation in X-rays, whilst the accretion disks in accreting 
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white dwarfs (cataclysmic variables or CVs) emit mostly 
in the optical/UV wavebands. This is a consequence of the 
gravitational potential well created by the central compact 
object: for CVs the inner-most edge of the accretion disk sits 
at a few thousand gravitational radii, whilst for galactic BH 
and NS the inner disk reaches down to a few gravitational 
radii. 

Most CVs and XRBs are highly variable sources in 
X-rays and/or optical/UV. This variability has been as- 
sociated with the accretion disk from the characteris- 
tic frequencies observed in the power spectral density 
(PSD). For example, the presence of periodic modula- 
tions slightly longer than the orbital period can be asso- 
ciated with positive superhu mps caused by a tidal defor- 
mation in the accretion disk (|WhitehurstJI 19881 : IOsakilll989l : 



2 S. Scaringi et al. 



lLubowlll99ll) This has been modelled and observed in CVs 



Montgomervll2012l; IWood et al-llioTll ), as well as in XRBs 
O'Donoehue fc Charleslll996l ). 



Another example of similar timing characteristics be- 
tween CVs and XRBs is the presence of both dwarf nova os- 
cillations (DNOs) and quasi-periodic oscillations (QPOs). In 
both cases DNOs/QPOs appear in the P SD at a few tens of 
mHz l|Warnerll2003 ; IPretorius et alj|2006h . whilst for XRBs 
they appear at few hectoHz, an d are referred to as lower kHz 
oscillations l|Belloni et al.1120021 ) . The phenomenological sim- 
ilarity between the QP Os and DNOs o bser ved in CVs an d 
XRBs was first noted bv lMauchd (|2002h and I Warner! l|2004 ). 
where the ratio of periods is Pqpo / Pdno ~ 15, and holds 
over 6 orders of magnitude in temporal frequency. The phys- 
ical reason for this relation is not fully understood, althoug h 
some suggestions have been proposed (|Belloni et al.l f2002). 
all involving accretion disk dynamics, but no ne seem to be 
able to explain all ob servations consistently (I Warner! |2004| ; 
iPretorius et~alll2006l ). 

In addition to periodic or quasi-periodic signals, both 
CVs and XRBs possess an intrinsic aperiodic broad-band 
noise continuum, generally described as flickering in CVs. 
This red noise component produces a continuum that rises 
towards low frequencies. In the XRB context, this is usu- 
ally modelled a by simple power-law (e.g. l/.f) or by a 
sum of broad Lorentzian components (|Belloni et al. 1 120021 1. 



with a characteristic break at 



10" J Hz for CVs and 



1 Hz for XRBs l|Revnivtsev et al.1 120121 ; iBelloni et all 
2002). This component has been shown to display a lin- 
ear rms-fl ux relation in XRBs and Active Galactic Nu- 
clei (AGN. IUttlev fc McHardvll200ll ), where the root-mean- 
square (rms) variability linearly scales with flux over a wide 
range of timescales. More recently, the same relation has 
been also observe d in the CV system MV Lyrae, albeit at 
lower frequencies (|Scaringi et al.l [20121 . hereafter Paper I). 
The detection of the rms-flux relation in both CVs and 
XRBs strongly suggests that the broad-band variability orig- 
inates within the accretion disk. Under certain assumptions 
about the origin of the variability and the physics of the 
mass-transfer process, the characteristics of the broad-band 
noise can be used to constrain the strength of the effec- 
tive viscosity that drives the mass transfer through the disk 
l|lngram fc Donell201Cl . l201ll ). 

Motivated by the similarities between the periodic and 
aperiodic signals in CVs and XRBs, we here study the broad- 
band temporal frequency spectrum of the CV MV Lyrae 
in search of further similarities between the two classes of 
objects. MV Lyrae is one of 14 known CVs in the Kepler 
field-of-view (FOV), and is classified as being a VY Scl no- 
valike system, spending most of its time in a high state 
(V ~ 12 — 13), but occasionally (every few years) under- 
going short-d uration (weeks to m onths) drops in brightness 
(V « 16 - 18, iHoard et al.1l2004l ). The reas on for these sud- 
den d rops in luminosity is not clear, but iLivio fc Pringld 
l| 19941 ) have suggested star spots from the donor covering 
the LI point inhibiting mass transfer. It is known however 
that MV Lyrae has an extremely low mas s transfer rate at its 
minimum br ightn ess (3 x 10~ 13 M^ Ivr, iLinnell et~aT]|2005l; 
IHoard et all 120041 '). where the WD is detected at V » 18 
and dominates the emitted light. Furthermore, the orbital 
period of 3.19 hours has been determined for the system, as 



well as a low inclination of i ~ 11° — 13° (Skillman et al. 
Il995h . 

In Paper I we analysed the high frequency (tens of min- 
utes) broad-band variability for MV Lyrae with data ob- 
tained with the Kepler satellite. In this paper we will anal- 
yse the same lightcurve at lower temporal frequencies, with 
particular emphasis on the varying broad-band components 
during the observation. MV Lyrae has been observed with 
Kepler during a low-to-high luminosity transition, when its 
optical emission originates mostly from its nearly face-on 
disk. 

In section O we briefly describe the Kepler data acqui- 
sition, and the procedure we use to construct broad-band 
power spectral densities (PSDs). We fit the PSDs with a 
combination of Lorentzian shaped functions to characterise 
the observed broad-band noise components. Section [3] shows 
our results, with particular emphasis on one frequency- 
varying broad-band noise component, which appears to cor- 
relate with the mean source flux. Finally, in section [4] we 
consider whether the inferred broad-band noise components 
can be associated with fluctuations on either the dynamical 
or viscous time-scales in the accretion disk. As we shall see, 
both of these interpretations lead to uncomfortable implica- 
tions for the accretion disk structure. 



2 DATA ANALYSIS 

The MV Lyrae lightcurv^B was provided to us by the Sci- 
ence Operations Centre in reduced and calibrated form af- 
ter being run throug h the standard data reduction pipeline 
(I Jenkins et al.l 12010). Similarly to Paper I, we only con- 
sider the Single Aperture Photometry (SAP) lightcurve. 
Data gaps occasionally occur due to Kepler entering anoma- 
lous safe modes (see the Kepler Data Characteristics Hand- 
bool0). Here, we make no attempt to correct these arti- 
facts, but simply remove them from the li ght curve. Fig. 
T] sho ws the short cadence (58.8 seconds; iGilliland et al.l 
2010J), barycentre corrected, lightcurve for MV Lyrae ob- 
tained during the first eight quarters of Kepler operations. 
The lightcurve spans an interval of 633 days. The visible 
data gaps present in the lightcurve are due to the artifacts 
described in Paper I as well as the monthly data down-links. 

In order to obtain the time-averaged PSDs, we have 
split the lightcurve in Fig. [T] into 120 non-overlapping seg- 
ments, each covering 5.275 days. This timescale has been 
chosen to probe the lowest frequencies without being af- 
fected by low frequency power generated by the long term 
trends in flux during the observation (see Fig. [TJ. We com- 
puted the Fast-Fourier transform (FFT) of 118 segments, 
ignoring 2 segment due to the large data gap at ~ 750 
d ays in Fig. [T] We furth er applied the rms normalisation 
of Miyamoto et alj |l99 ll ) so that the square root of the in- 
tegrated PSD power over a specific frequency range yields 
the root-mean-square (rms) variability. After computing the 
normalised PSDs for each segment we obtain the intrinsic 
spread in each frequency bin by computing the standard 
deviation of the 118 PSDs, shown in Fig. [2] 



1 The Kepler response function for the photometry covers the 
wavelength range 4000 - 9000A. 

2 http:/ /archive. stsci.edu/kepler/documents. html 
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Figure 1. MV Lyrae light curve obtained by the Kepler satellite in short cadence mode (58.8 seconds). 
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Figure 2. Time-averaged PSD of MV Lyrae. The shaded region 
shows the 1-ct width of the distribution spanned by the 118 single 
5.275 day PSDs, whilst the thick line shows their average. 
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Figure 3. Rms-normalised dynamic PSD of the MV Lyrae 
lightcurve shown in Fig. [T] Each PSD has been computed on 
a 5.275 days segment, with a 50% overlap. 



At first glance we note large scatter in the time-averaged 
PSD, caused by the varying PSD shape and normalisation 
during the 633 day observation. Furthermore many incoher- 
ent features are clear from Fig. [5] Specifically, we see a clear 
broad Lorentzian-shaped feature peaking at about 10 -3 Hz, 
briefly discussed in Paper I, but also excess power at about 
1CP 4 and 1CF 5 Hz. It is hard to determine whether the ex- 
cess power is caused by a single intrinsically broad feature or 
a blend of many. In order to shed light on which of these in- 
terpretations is correct we have produced a dynamic FFT of 
the lightcurve. This is shown in Fig. [3] where we have again 
computed rms-normalised PSDs on 5.275 day segments, but 
this time with a 50% overlap, for clarity. The apparent gaps 
are caused by data gaps in the lightcurve shown in Fig. 
H Fig- E] suggest that the excess power between 10 -3 ' 5 Hz 



and 10 -4 Hz is produced by a single, frequency-varying QPO 
whose peak frequency is initially near 10 _3,5 Hz and then 
slowly declines on a time-scale of » 250 days. A comparison 
with the lightcurve in Fig. [T] also suggests a possible cor- 
relation between the peak frequency of this QPO and the 
mean flux of the system. We will return to this correlation 
in section [3] 

2.1 PSD fitting 

In order to quantify and fit the frequency- varying QPO(s), 
we first produced 118 rms-normalised PSDs, each from 
5.275-day, non-overlapping data segments. Next, we aver- 
aged 5 consecutive PSDs (except in one case where we aver- 
aged 3 due to a large data gap in the lightcurve) and binned 
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these into ~ 550 equally spaced frequency intervals. The 
bin size is somewhat arbitrary, but has been chosen as a 
compromise to increase the signal-to-noise in each bin while 
retaining a high frequency resolution. We further obtained 
the intrinsic scatter in each frequency bin by calculating the 
standard deviation of the 5 PSDs used. This intrinsic scatter 
is larger than the statistical errors, and we thus adopt the 
intrinsic scatter as a conservative estimate of the errors of 
each frequency bin, which also takes into account fluctua- 
tions in the intrinsic PSD shape. This procedure resulted in 
24 time-averaged PSDs. 

The fit to each PSD was carried out by starting with a 
simple broken power-law fit, to which we added successive 
Lorentzian component until the reduced \ 2 w as near unity. 
The s hape of the Lorentzians was taken from Belloni et al.l 



P(u) 



TV A 2 + (V-V ) 



(1) 



where r is the integrated fractional rms and A is the half- 
width at half-maximum. With this definition the Lorentzians 
attain their maximum power in vP(y) at the characteristic 
frequency 



^ 2 +A 2 . 



The broken power-law takes the form 



1 + 



(2) 



(3) 



where A is the normalisation and ut r the break frequency. 
With this definition the PSD shape has a slope of —1 up 
to ifrr, where the slope then changes to —2, similarly to the 
PSD s observed in some AGN and XRBs l)Summons et al.l 
120071 ). In our analysis we are mainly interested in the prop- 
erties of the fitted Lorentzian, and use the low-frequency 
power-law model to account for any systematic excess power 
generated by the long term trends (such as the steady flux 
increase) observed in the lightcurve shown in Fig. [T] 
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Figure 6. Correlation between the characteristic timescale t max 
of the Lorentzian L3 and mean source flux. 
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Figure 7. Evolution of the integrated fractional rms n of the 
Lorentzian L\ as a function of observing time. 



3 RESULTS 

We fit each of the 24 PSDs individually, taking into account 
the Poisson noise level. We find that all PSDs require 4 
Lorentzians. Of course, this does not exclude the possibility 
that more components could be present in the data, hidden 
beneath the noise level. Fig. [4] shows some examples of our 
fits, whilst Table [T] shows the resulting best-fit parameters 
for all our fits, where we labelled our Lorentzians from 1 to 

4 in ascending order of frequency. In a few cases when two 
Lorentzians were highly blended the fits resulted in some 
parameters being unconstrained. 

We searched for possible correlations within the fit 
parameters and both observation time and mean source 
flux. In this respect, interesting trends with observation 
time are found within the characteristic frequencies of the 
Lorentzians v max , shown in Fig. [5] The most prominent 
correlation however is found between the characteristic fre- 
quency of Lorentzian L3 and the mean source flux, which 
obtained a Spearman correlation coefficient of —0.89, and is 



shown in Fig. [6] To date, no similar correlation to that ob- 
served here has been reported for a CV, and we will discuss 
the possible origin of such variability in the next section. No 
other fit parameters display a clear trend with either flux or 
time. 

The high-frequency break, modelled by the Lorentzian 
Li, is seen to slightly change in frequency, but does not 
correlate with mean flux (see top-left panel of Fig. [5} • How- 
ever, one interesting trend within the Lorentzian L\ arises 
by looking at the integrated fractional rms, ta, as a func- 
tion of time, shown in Fig. [7] The two peaks at « 250 and 
~ 500 days coincide with a small plateau and a decrease 
in flux within the MV Lyrae lightcurve respectively. Both 
these features, and the break frequency, will be discussed in 
more detail in the next section. 

As mentioned in the introduction, correlations between 
some of the Lorentzian components are observed in XRBs 
and some CVs. Specificall y, the two-QPO co rrelation dia- 
gram produced for XRBs ((Belloni et alj|2002h and for CVs 
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Figure 4. PSDs in vP v form for MV Lyrae, also showing the best-fit model components. The PSDs have been created by averaging 5 
5.275 day PSDs. From top-left to bottom-right the central observation times are respectively 182.7, 367.3, 525.5, 789.1 in BJD-2454833. 



\ Warner! [2004) can be described by Pqpo/Pdno ~ 15. We 
have searched for a similar correlation between the 4 fitted 
Lorentzians in MV Lyrae. In Fig. [8] we show plots for the 
characteristic frequencies of all four fitted Lorentzians. In 
order to determine whether any of these pairs is statistically 
correlated we perform ed a Monte-Car l o corr elation analysis 
similar to that used in IScaringi et all j2QlcT >. 

For each pair of characteristic frequencies we created 
two datasets. In the first (set A) we replaced each of the 
values with a random variable drawn from a normal dis- 
tribution whose mean is equal to the fitted characteristic 
frequencies and whose standard deviation is equal to the re- 
spective errors. In the second (set B) we additionally shuffled 
the order of one of the two frequencies, thus randomising 
the correlation rank. We produced 200,000 samples for each 
set, and computed the correlation coefficient p for each. The 
mean values of the p distributions obtained from set A will 
then yield the intrinsic correlation coefficients, taking into 
account the errors. The corresponding set B distribution will 
instead yield the width of the p distribution in the absence 
of any correlation (which is centred at 0). The statistical sig- 
nificance for each pair is then found by determining where 
the obtained mean p values from set A fall on the obtained 



distribution of our mock dataset B. This analysis thus allows 
us to determine the correlation significance of the observed 
characteristic frequency pairs by taking into account the er- 
rors on each observation. 



In Fig. [S] we show in the top-right panels plots for the 
characteristic frequency pairs, whilst in the bottom-left we 
show the resulting p distribution from our first simulation 
and the corresponding correlation significance. None of the 
Lorentzian pairs are observed to display a clearly statisti- 
cally significant correlation. However, it is still possible that 
some pairs may be correlated, and the low-frequency clus- 
ter of datapoints is biasing our observation. For example, 
the pair v max ^ vs. Vmax,2, which obtained the highest sig- 
nificance (<r = 2.6) in our analysis, as well as some other 
pairs, do suggest this. However, further measurements are 
required at the highest frequencies in order to determine if 
this is the case. If the v ma x,3 vs. v max $, pair were to be 
correlated after all, our analysis suggests they would follow 
the relation t max ,3/tmax,2 ~ 3, different from the known 
Pqpo /Pdno ~ 15 relation (|Belloni et all |2002| ; IWarnerl 
l2004h . 
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20 


9 


02+0 


61 





17+0 


08 


5.55+0.44 


9.89+0.51 


0.10 


491 


48(540) 


2 


95+0.20 


1 


47+0 


01 





87+0 


11 


8 


89+2 


06 





24+0 


16 


4.60+0.47 


5.39+37.38 


0.18+1.59 


620 


69(539) 


2 


78+0.22 


1 


39+0 


01 


1 


10+0 


31 


8 


41+0 


38 





40+0 


03 


5.67+0.54 


15.15+0.48 


0.10 


542 


88(540) 


2 


35+0.13 


1 


27+0 


02 





99+0 


14 


2 


76+0 


19 





59+0 


01 


1.35+0.31 


1.47+0.03 


68.14+11.31 


646 


78(539) 


1 


68+1.90 


1 


26+0 


05 





29+0 


99 


7 


48+0 


76 





40+0 


08 


5.52+0.52 


7.02+2.73 


0.46+0.31 


652 


78(539) 


3 


00+0.27 


1 


24+0 


01 


1 


14+0 


30 


6 


34+0 


37 





57+0 


02 


3.45+0.44 


3.80 


1.18+0.09 


571 


55(540) 


3 


84+3.33 


1 


22+0 


01 





34+0 


89 


8 


44+1 


31 





43+0 


07 


4.07+0.66 


5.01+1.35 


1.14+0.92 


611 


36(539) 


2 


48+0.47 


1 


46+0 


04 


1 


88+0 


63 


7 


68+0 


36 





55+0 


01 


1.87+0.15 


11.23+1.64 


0.72+0.21 


503 


37(539) 


3 


17+0.41 


1 


58+0 


04 


2 


73+0 


77 


9 


50+0 


59 





45+0 


02 


2.51+0.25 


15.80+4.22 


0.48+0.22 


654 


27(539) 


1 


99+0.63 


1 


87+0 


08 


2 


31+1 


45 


7 


78+0 


93 





60+0 


08 


5.35+0.96 


7.09+1.44 


0.82+0.37 


663 


23(539) 


2 


41+0.24 


1 


60+0 


03 


2 


34+0 


41 


4 


62+0 


27 





64+0 


02 


2.98+0.35 


3.41+0.40 


0.90+0.23 


561 


41(539) 


3 


27+0.53 


1 


34+0 


04 


2 


85+0 


63 


5 


70+0 


89 





54+0 


05 


3.75+0.95 


1.46+1.24 


0.57+0.94 


657 


11(539) 


1 


92+0.23 


1 


42+0 


02 


1 


21+0 


30 


6 


68+0 


32 





55+0 


03 


5.31+0.50 


1.90 


0.47+0.03 


459 


03(540) 


2 


40+10.08 


1 


39+0 


02 





21+2 


01 


9 


51+1 


15 





24+0 


12 


6.18+0.47 


1.02+6.46 


0.32+3.10 


589 


60(539) 


1 


90+0.59 


1 


32+0 


04 


1 


94+0 


94 


7 


52+0 


81 





48+0 


06 


5.35+0.80 


1.50+0.66 


0.96+1.12 


556 


03(539) 


2 


23+3.48 


1 


38+0 


01 





30+1 


34 


5 


51+0 


57 





75+0 


03 


3.32+0.56 


1.06+0.09 


5.48+4.20 


565 


27(539) 


3 


20+13.70 


1 


54+0 


01 





15+1 


38 


8 


13+0 


89 





36+0 


07 


4.22+0.36 


5.01+16.93 


0.19+0.87 


602 


86(539) 


2 


21+0.34 


1 


41+0 


03 


1 


74+0 


55 


8 


29+1 


13 





32+0 


11 


5.30+0.54 


5.01+17.83 


0.19+0.89 


538 


07(539) 


3 


21+0.63 


1 


34+0 


04 


2 


29+0 


70 


5 


77+0 


72 





68+0 


03 


3.55+0.90 


10.56+6.00 


0.21+0.16 


562 


37(539) 



Table 1. Best-fit parameters for the 24 time-averaged PSDs of MV Lyrae. In cases where the fit parameters were unconstrained no error 
bars are given. 



4 DISCUSSION 

Here we discuss the implications of the observed 
QPO/broad-band noise component frequencies with regard 
to the stan dard, optically thick, geom etrically thin accretion 
disk model (jShakura &i Sunvaevll973p. Given the low binary 
inclination of i = 11° - 13° (jSkillman et al.lll995l ). the WD 



magnitude seen in the faint state l|Hoard et al.l [20041 ) and 
the result from Paper I, we are confident that most of the 
emitted light observed in the Kepler passband is originat- 
ing from the accretion disk. In what follows w e will always 
assume M WD = 0.722M Q (|Hoard et al.ll2004l ), and a WD 
radius of w 7 x 10 8 cm. 

The observed characteristic frequencies are usu- 
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Figure 5. Characteristic frequency, Umax, versus observation time for the four Lorentzians L\, L2, L3 and L4 from top-left to bottom- 
right respectively. 



ally associated w it h the dynamical or vis c ous t im escale 
llRevnivtsev et al.1 120121 : llngram fc Pond 120101 , l201ll ; 
lUttlev fc McHardvll200lf '). While it is not yet clear if one 
of these two timescales is the correct one, most other 
time-scales of interest (such as thermal) lie between these 
two. We will thus discuss them as likely limiting cases. The 
dynamical timescale is given by 



GM 



tdyn{r) V r 3 4n 2 ' 
and relates to the viscous timescale, 

Vvisc{r) = a(H/R) 2 u dyn {r), 



(4) 



(5) 



where r is the accretion disk radius, H/R and a are the disk 
scale height and viscosity parameter respective ly, and M is 
the m ass of the central compact object (see e.g. iFrank et al.l 
l2002h . 

For XRBs and AGN the observed broad-band 
noise variability is generally asso ciated with the viscous 
timescales over a range of radii (|Arevalo fc Uttlevl [2006; 
llngram fc Pond l20ll l201lf ). This is motivated by obser- 
vation s of the rms-flux relation seen in many XRBs and 
AGN (|Uttlev fc McHardvll200ll ). which imply that the fast 



variability must be coupled to the slow variability by multi- 
plicative processes (i.e. the resulting variability at the inner- 
most edges of the disk is the product of the variability pro- 
duced at larger annuli). This can naturally arise as a re- 
sult of the viscous propagation of material inwards through 
the disk. The standard model which seeks to reproduce 
these observations is the fluctuating accretion disk model 
l|Kotov et al.l l200ll ; iLvubarskiil Il997l ) , which can correctly 
reproduce the rms-flux relation as well as the general 1// 
shape and high-frequency break observed in most PSDs of 
XRBs and AGN. In both cases the resulting inner-edge of 
the disk is inferred as being geometrically thick in order 
to allow the variability to propagate inwards without be- 
ing damped. On the other hand, some authors have instead 
associated the break fre quency observed in CVs to t he dy- 
namical timescale (Eq|4l [Revnivtsev et al.|[2010l . l2012h of the 
inner-most edge of the disk. This association seems to alle- 
viate the problem of a thick disk (since for this interpreta- 
tion the disk can be either thin or thick), but instead sug- 
gests a l arge disk truncation radius at IQ Rwd , and evap- 
oration l|Mever fc Mever-Hofmeisterl 1 1994 ) is invoked as a 
mechanism to deplete the inner disk. Recentl y the rms-flux 
relation has also been observed MV Lyrae ( Scaringi et al.l 
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Figure 8. The top-right half of the figure shows the characteristic frequencies of the fitted Lorentzians. The bottom-left part of the 
figure shows the corresponding Spearman rank correlation coefficient values p obtained from the analysis described in section [3] together 
with the corresponding correlation significance. Each panel matches its reflection along the diagonal. 



2012) with the same data as that used in this work, suggest- 
ing that the mechanism responsible for mass transfer within 
the disks, and the associated variability, is similar for CVs, 
XRBs and AGN, whether this is viscous, dynamical or ther- 
mal. The main difference between the rms-flux relation ob- 
served in XRBs and that observed in the CV MV Lyrae are 
the timescales involved: for XRBs the rms-flux relation is 
observed in X-rays on timescales of seconds whilst for MV 
Lyrae it is observed in optical light on minute timescales. 
This timescale difference can potentially be attributed to 
the different emitting regions, since the variability observed 
in X-rays in XRBs is originating deeper within the potential 
well (thus faster variability) of the accreting compact source 
when compared to MV Lyrae. The study performed in Pa- 
per I concentrated on analysing the high-frequency proper- 
ties of MV Lyrae, whilst here we concentrate on the broad- 



band variability properties, and assess whether the observed 
features are phenomenologically similar to those in XRBs. 

Given that the nature of variability is still debated, but 
recent phenomenological similarities seem to suggest a com- 
mon origin for CVs, XRBs and AGN, we will, in what fol- 
lows, associate the characteristic frequencies observed in MV 
Lyrae to either the viscous or dynamical timescales, and in- 
fer values for the accretion disk structure. This will allow us 
to compare our results to those inferred from X-ray studies 
of XRBs and AGN. In the next subsection we will specifi- 
cally concentrate on the dynamical/ viscous interpretations 
of two Lorentzians, namely Li and L3. The general conclu- 
sions drawn from these two components will also hold for 
the other two Lorentzians as well, and we choose to discuss 
only L\ and L3 as these components seem to display the 
most interesting trends in Fig. [5] [6] and [7] 
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4.1 Viscous interpretation 



4- 1.1 The low frequency Lorentzian L3 

Here we discuss the possibility that the observed characteris- 
tic frequency v m ax,3 is associated with the viscous timescale 
(t v i sc ) at some disk radius. The apparent anti-correlation 
with flux seen in Fig. [6] could then be suggesting a change 
in the viscosity parameter a and/or scale height H/R. Ob- 
servationally, we find that the low frequency Lorentzian L3 
produces an rms variability of ss 5%. As at least that frac- 
tion of the luminosity needs to be emitted in the variable 
emission region, we can first obtain a lower limit on the size 
of this emitting region. 

Measured temperature profiles of CVs with a similar 
accretion rate seem to follow the theore tical predictions o f 
standard accretion discs up to w 0.057?© (|Groot et alj|200ll) . 
Closer in, there is some controversy on the shape of the tem- 
perature profile, especially for the subcl ass of CVs SW Sex 
stars of which MV Lyrae is a member of. iGroot et al.l |2004l . 
l200ll ), have determined a flat-top profi le thro ugh in direct 
eclipse mapping, however iKnigge et alj l|2000l . [2004) have 
later confirmed the possibility that in the high-luminosity 
state, at least some SW Sex stars posses a self-occulting 
component. This could be due to the presence of a thick 
disk, and would thus result in an observational bias when 
measuring the temperature profiles in the inner-most regions 
of the accretion disk, producing a flat-top temperature pro- 
file. 

Because of these controversies we will take the conser- 
vative limit and model the minimal size of the emitting re- 
gion by using the standard disk temperature profile from 
iFrank et all i|2002h ( Eg. 5.49) wit h a mass accretion rate 
of 3 x lO" 9 M 2/r -1 ijLinnell et al.ll2005h . We then summed 
the emitted blackbody spectra produced at each disk an- 
nulus and folded this with the Kepler response functior0, 
producing the relative cumulative flux emitted by the disk 
as a function of radius in the top panel of Fig. [9] The fig- 
ure shows, for example, that if we observe say « 5% rms 
variability (as observed for L3), the emitting region must 
be larger than « 3Rwd- In fact, 3Rwd in this scenario is a 
conservative lower limit for two reasons: (i) the absolute am- 
plitude of the variabil ity is larger than s a 5%, and (ii) some 
CVs (namely SW Sex. lGroot et al.ll200ll ) show that the disk 
temperature profile flattens out at larger radii (~ O.2.R0). 
Both these effects will push the minimal size of the emitting 
region to larger radii than those shown in Fig. [9] and further 
increase the resulting limits on a and H/R. 

Having produced an estimate for the minimal size of 
the variable emitting region, we now place limits on the 
disk scale height (H/R) and viscosity (a) for the fastest 
measured timescales of the Lorentzian L3 of 0.8 hours. This 
is done as a function of radius in the bottom panel of Fig. 
[9] (solid lines). The inferred values for a and H/R axe too 
high to be accommodated by the standard thin disk model, 
where a ^ 0.1 and H/R ^ 10~ 2 are usually taken. Addi- 
tionally, the most conservative values for a and H/R (solid 
black line) are ruled out by the constraints imposed by the 
minimal size of the emitting region. Thus, if one wants to 
adopt the viscous interpretation, our result seems to sug- 
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Figure 9. Top: Cumulative radial emission profile for the ac- 
cretion disk in MV Lyrae folded through the Kepler passband. 
Bottom: Inferred accretion disk radius and viscosity parameter 
a for different disk scale heights. The estimates have been ob- 
tained by associating the characteristic frequencies of 0.8 hours 
(L3, solid lines) and 13 minutes (Li, dashed lines) to the viscous 
timescale at a particular disk radius. The inferred a and H/R val- 
ues are found to be too large to be accommodated by a standard 
thin disk. The vertical black dashed lines indicate the minimal 
extent of the accretion disk given the observed 5% and 10% rms 
amplitude of the variability. 



gest the presence of a geometrically thick accretion disk in 
order to explain the observed characteristic frequencies of 
the L3 Lorentzian. This wou ld be in agreement with the 
work o f[Ch urazov et al.ll|200j ). who have noted that for the 
fluctuating accret ion disk model to work (|Lvubarskiilll997l ; 
iKotov et al.l 200"ll) and produce the observed rm s-flux rela- 
tion (|Arevalo fc Uttley|2006l : IScaringi et al.l l20I2l). one needs 
a geometrically thick disc. 



4-1.2 The high frequency Lorentzian L\ 

We find the characteristic frequency of the Lorentzian L\ 
to slightly vary throughout the observation, but the data 
does not allow us to determine whether it correlates with 
flux. In Paper I (Fig. 7) we inferred values of H/R « 0.1 for 
a ~ 1 assuming the ss 13 minute break originates from very 
close to the WD surface. Here we extend that analysis, and 
show in Fig. [9] (dashed lines) the dependence on both H/R 
and a for various disk radii using the inferred PSD break at 
about 13 minutes to represent the viscous timescale. In this 
case the minimal size of the emitting region is determined 
by the w 10% rms variability observed in Fig. Like the 
situation with the 0.8 hour Lorentzian (see section [4.1.11) . 
this last constraint rules out the most conservative values 
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for a and H/R, suggesting a geometrically thick disk. Fur- 
thermore we note that the constraints inferred from Fig. [9] 
are likely to be even larger than those inferred here. From 
the contin uum variability modelling of accretion disks per- 
formed by llngram fc Pond (|201Cj| . l201ll ) we know that the 
true viscous timescale at the inner-most edge of the disk is 
even faster than that observed from the PSD break, pushing 
the inferred a and H/R to even higher values. 

One other important result concerning the high- 
frequency Lorentzian L\ is the change in fractional 
rms variability shown i n Fi g. [Jj From the work of 
iGeertsema fc Achterbergl (|l992l ). we can use the fractional 
rms as an indication of both a and H/R. This method is 
based on the effects of magneto-hydrodynamic (MHD) tur- 
bulence in an accretion disk, and associates the source of 
flickering to the energy dissip ated per unit area at t h e disk 
surface. With this method iBaptista fc Bortolettol (|2004l ) 
have inferred a « 0.16 for H/R = 1CP 2 in the CV sys- 
tem V2051 Oph, and comment on this result as being un- 
comfortably high to be accommodated by the standard thin 
disk model. In Paper I we also showed how a varies between 
w 0.16 and ~ 0.62, and invoked an even thicker disk of 
H/R ~ 0.1 to avoid at least the highest a values. The peaks 
of the rms variability presented in Fig. [JJ coincide with a 
short-duration (few weeks) plateau observed at « 250 days, 
and a drop in flux seen at ~ 500, both shown in Fig. [T]). The 
behaviour observed in Fig. [JJcould be suggesting a change in 
the geometry and/or viscosity parameter during the obser- 
vation since t he rms variability is tied to both a and H/R 
in the work of IGeertsema fc Achterberd (| 19921 ). If this were 
the case however, the inferred range for both a and H/R is 
quite large (0.02 < a < 0.2 and 0.01 < H/R < 0.1) during 
the rms change between 5% and 10%, and we caution inter- 
pretation of this result as it is not yet clear what physical 
mechanism is required to explain the observed phenomena. 



4.2 Dynamical interpretation 

4-2.1 The low frequency Lorentzian L3 

We can can alternatively try to associate the obtained char- 
acteristic frequencies of the Lorentzian L3 from Fig. [5] with 
the Keplerian frequency at a specific disk radius. In order 
to explain the correlation seen in Fig. [6l the emitting region 
would then radially grow in size, and increase in tempera- 
ture, as the timescale changes from 0.8 to 2.1 hours. 

Using the highest and lowest frequencies measured for 
the Lorentzian L3, we can infer using Kepler's third law an 
accretion disk size ranging from ~ 0.39i?o, for the highest 
frequency, up to ~ 0.7ARq for the lowest. These estimates, 
especially the low frequency one, are not consistent with the 
estimated position of the L I point in this bi nary system, 
which sits at about 0.577?© (|Frank et alj|2002h . One way to 
overcome this inconsistency would be to assume that the ob- 
served frequency is a beat between the Keplerian outer disk 
frequency (fdisk) and the orbital frequency (p or b) of 3.19 
hours, Vbeat = i^disk + Vorb- ln this case, the outer disk-edge 
would be observed to grow from « O.33i?0 to « O.53i?0, 
consistent with the disk radially growing in size up to the 
LI point. However the other Lorentzians, especially the low- 
est frequency L4 Lorentzian, would then require a different 



mechanism to drive them, as their beat frequency would be 
too low to be accommodated within the LI point. 



4-2.2 The high frequency Lorentzian L\ 

The high-frequency PSD breaks observed in accreting com- 
pact objects have also been associated with the Keplerian 
timesc ale at the inner-most ed g e of a ccretion disks. Specif- 
ically, iRevnivtsev et all (|2012l , |2010| ) have associated the 
PSD breaks observed in some magnetised WDs, NSs and 
the well-known CV SS Cyg to the dynamical timescale 
at the inner-most edge of the accretion disk. They in- 
ferred relatively large inner disk truncation radii for all sys- 
tems they studied, and for SS Cyg in particular they in- 
ferred a truncation radius of fa lORwp ■ In order to deplete 

suggest evaporation 



the inner disk Revnivts ev et al.l .._ , 
|Mever fc Mever-Hofmeistei 19941 ^1 as a mechanism to trun- 



cate the thin disk away from the WD surface. 

The break frequency observed in SS Cyg is very similar 
to the one observed here in MV Lyrae, and it is possible that 
we also observe the Keplerian frequency at the inner-most 
edge of the accretion disk truncated at about ~ IORwd- We 
note, however, that dynamical processes will be more easily 
damped, as compared to viscous process, whilst they prop- 
agate through disk. This is because viscous damping will 
suppress the fast dynamical variability but not the slower 
viscous variability. As a result dynamical effects are observed 
as additive process in lightcurves, whilst viscous ones as mul- 
tiplicative. The presence of the rms-flux relation at high 
frequencies thus suggests viscous interactions as the main 
source of variability, since additive processes would destroy 
the observed rms-flux relation. 



5 CONCLUSION 

We have presented an analysis of the broad-band frequency 
behaviour of the accreting WD MV Lyrae based on data ob- 
tained by the Kepler satellite. We have shown how the com- 
plex PSD can be decomposed with a number of Lorentzian- 
shaped functions. We further searched for possible correla- 
tions between the characteristic frequencies, and found the 
first frequency varying QPO in a CV, where frequency is 
inversely proportional to mean source flux. 

The characteristic frequencies associated with the fit- 
ted Lorentzians were used to explore the origin of variable 
emission in terms of viscous or dynamical processes as the 
two limiting cases. In the former case we infer extremely 
high values of both disk scale height H/R 0.3, and vis- 
cosity a ^ 0.1, suggesting the existence a geometrically- 
thick disk. This result is p otentially in line with the work of 
iKnigge et all |2000l . 120041 ), which have suggested the pres- 
ence of a self-occulting accretion disk close to the WD. In 
the dynamical case we instead infer a large disk truncation 
radius of ~ IORwd, but the presence of other components 
undermines a dynamical interpretation, at least for the low- 
est frequency components. More importantly, the presence of 
the r ms-flux relation observe d in MV Lyrae ([Scaringi et al.l 
I2012D and other XRBs/AGN (|Uttlev fc McHardvll200ll ). po- 
tentially rule out dynamical effects as the source of variabil- 
ity (since the fast dynamical variability will be damped), but 
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favours a viscous origin for the observed broad-band noise 
components, at least at the highest observed frequencies. 

In summary both viscous and dynamical (and conse- 
quently thermal) timescales struggle to consistently explain 
the observed broad-band variability in MV Lyrae. However, 
our analysis seems to suggest, at least phenomenologically, 
that the mechanisms which give rise to the observed rms- 
flux relation(s) and characteristic frequencies observed here 
and in other XRBs/AGN need to occur in all accretion discs 
(whether they are thin or thick), possibly suggesting a sim- 
ilar physical origin for the variability in both types of sys- 
tems. 

The broad-band variability properties of MV Lyrae re- 
main yet to be fully understood. More generally, the vari- 
ability properties of both CVs and XRBs also remain an 
enigmatic observational feature. Although both types of sys- 
tems possess an accretion disk which is in many ways similar, 
there has not yet been enough observational data to provide 
a complete comparative study of the broad-band variability 
properties. However, it is clear that the phenomenological 
properties between the broad-band variability observed in 
MV Lyrae and in X-ray binaries and Active Galactic Nu- 
clei are very similar (when appropriately scaled), and are 
potentially driven by a common accretion mechanism. 
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